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Abstract 

Let n be an open period annulus of a plane analytic vector field 
Xq. We prove that the maximal number of limit cycles which bifurcate 
from n under a given multi-parameter analytic deformation X\ of Xq 
is the same as in an appropriate one-parameter analytic deformation 
^A(e)i provided that this cyclicity is finite. Along the same lines we 
give also a bound of the cyclicity of homoclinic saddle loops. 

1 Statement of the result 

Let Xx, A G (M", 0) be an analytic family of plane vector fields, and let 11 be 
an open period annulus of Xq (an open domain which is a union of periodic 
orbits of Xq). The cyclicity Cyclijl, X\) of 11 with respect to the deformation 
Xx is the maximal number of limit cycles of Xx which tend to 11 as A tends 
to zero, see Definition [3] bellow. The reader should not confuse the cyclicity 
Cyclijl, Xx) of the open period annulus 11 with the cyclicity Cycliti, Xx) of 
the closed period annulus fl. 

Our first result is the following 

Theorem 1 // the cyclicity CycKjl^Xx) of the open period annulus 11 is 
finite, then there exists a germ of analytic curve e ^ A(£), e G (M, 0), A(0) = 
0, such that 

Cycl{Il,Xx) = Cycl{Il,Xxi,))). (1) 
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In other words, the problem of finding the cyclicity of an open period annu- 
lus with resect to a multi-parameter deformation, can be always reduced to 
the "simpler" problem of finding cyclicity with respect to a one-parameter 
deformation. Indeed, in this case the displacement map can be expanded 
into power series 

d{u,e) =e''{M^{u) + eR{u,e)) (2) 

where Mg, the so called higher order P oincare-P ontryagin (or Melnikov) func- 
tion, depends on the germ of analytic curve ^ : e ^i^)- If A = (0, 1) is an 
interval parameterizing a cross section to the annulus 11, then is analytic 
on (0, 1) and its number of zeros Z{M^) counted with multiplicity is an upper 
bound of Cyd(Jl,Xx{e)))- Therefore we get the inequality 

Cyd{U, Xx) < sup Z{M^) (3) 

where the sup is taken along all germs of analytic curves 

e:£^A(£),e(0) = 0. 

Let X\ be an arbitrary analytic deformation of a vector field Xq. It 
has been conjectured by Roussarie [20l p. 23] that the cyclicity Cyd{T, Xx) 
of every compact invariant set T of Xq is finite. A particular case of this 
Conjecture is therefore that Cydijl, Xx) < oo, that is to say the claim of 
Theorem [T] holds without the finite cyclicity assumption on the open period 
annulus. It follows from the proof of Theorem [T] that if Cyd(Jl, Xx) = oo, 
then there exists a germ of analytic curve e i— >• X{e),e G (M, 0), A(0) = 0, 
such that the corresponding higher order Poincare-Pontryagin function M^{u) 
defined by ((21) has an infinite number of zeros in the interval A. In the 
particular case when Xq has an analytic first integral in a neighborhood of 
the closed period annulus IT, the analytic properties of M^{u) are studied in 
[HI [in]. Using this, it might be shown that M^{u) has a finite number of zeros 
on the open interval A. 

Suppose that the open period annulus of Xq contains in its closure a 
non-degenerate center. Denote the union of 11 with such a center by 11. If 
we blow up the center, it becomes a periodic orbit of a new vector filed to 
which Theorem [T] applies with minor modifications. It follows that ^ holds 
true with 11 replaced by 11. The question whether we can replace the open 
period annulus 11 by its closure 11 is much more delicate. Namely, suppose 
that the closed period annulus 11 is a union of 11, a non-degenerate center, 
and a homoclinic saddle connection (as for instance the two bounded annuli 
on figH]). By a homoclinic saddle connection (or separatrix loop) we mean 
a union of hyperbolic saddle point with its stable and unstable separatrices 
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which coincide. The union of these separatrices is a homoclinic orbit of the 
vector field, and together with the saddle point they form a separatrix loop 
or a homoclinic saddle connection. 

Suppose that A = [0, 1] parameterizes a cross section to 11. Consider 
as above a germ of analytic curve ^ and the corresponding higher Poincare- 
Pontryagin function M^. It is continuous on [0, 1], analytic on [0, 1) and has 
an asymptotic Dulac series at m = 1 

oo 

M^{u) = Y,(^2^{l-uy + a2i+i{l-uY+Hog{l-u) . (4) 

i=0 

We define a generalized multiplicity of a zero of at m = 1 to be equal to 
j if ao = fli = ■ ■ ■ = ttj-i = 0, aj ^ 0. At m = we define the multiplicity of 
the zero of 

oo 
i=l 

to be equal to j, where ai = ■ ■ ■ = = 0, a^+i 7^ 0. Define finally Z{M^) to 
be the number of the zeros (counted with multiplicity) of on [0, 1]. It is 
classically known [22j that 

Cyc/(n,XA(.)))<z(M5). 

We shall prove 

Theorem 2 Under the above conditions 

Cyd{Ti, Xx) < sup Z{M^) (5) 

where the upper bound is taken along all germs of analytic curves 

e:£^A(£),e(o) = o. 

Remarque. The hyperbolicity of the saddle point is essential for the proof 
of Theorem\^ It follows from [2l[ Proposition A 2.1, p. Ill] that the vector 
field Xq possesses an analytic first integral in a neighborhood o/II. Therefore 
Xq is a Hamiltonian vector field with respect to a suitable area form. It is 
also known that the cyclicity of the closed period annulus in this case is finite 
12^. 

Example. TheoremU\ applies to the three open annuli on figUi while Theo- 
remlE applies only to the two closed and bounded period annuli on figUi 

The proofs of Theorems [T] and [2] are inspired by the Roussarie's paper 
[T9] . in which the inequality ([5]) is shown to be true for a single regular peri- 
odic orbit. The equality ^ for a single regular periodic orbit is announced 
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Figure 1: Period annuli 



by Cauberg [3j (see also [4] for related results) . The finiteness of the cyclicity 
Cyd{'y, Xx) of a regular periodic orbit 7 was previously proved by Frangoise 
and Pugh [7]. The main technical tools in the proof of Roussarie's theo- 
rem [19| are the Hironaka's desingularization theorem applied to the Bautin 
ideal, followed by a derivation-division algorithm. This second argument ap- 
plies only locally. To prove Theorem [T] we also use the Hironaka's theorem, 
but replace the derivation-division algorithm by a variant of the Weierstrass 
preparation theorem. This already gives an upper bound of the cyclicity in 
terms of zeros of Poincare-Pontryagin functions. To get the exact result ^ 
we use the curve selection Lemma, as suggested by Roussarie [19]. Theorem 
[2] has a similar proof but is based on [22l Theorem C]. We note that we 
obtain a non-necessarily exact upper bound for the cyclicity of the closed 
period annulus. The reason is that the bifurcation diagram of limit cycles 
near the separatrix loop is not analytic, and therefore we can not apply the 
curve selection lemma. 

We mention finally that Theorem [T] end [2] allow an obvious complex ver- 
sion (with the same proof). For this reason most of the results in the next 
section are stated in a complex domain as well. The question about finding 
an explicit upper bound of the cyclicity of more complicated separatrix con- 
nections is almost completely opened. A recent progress in this direction is 
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obtained in [6j. 

The paper is organized as follows. In the next section we formulate several 
classical results which will be used in the proof. The latter is given in section 
[31 In the last section [4] we discuss some open questions. 

2 Digression 

In this section we formulate, for convenience of the reader, several facts of 
general interest, which are necessary for the proof of Theorem [U The base 
field is K = M or C. The corresponding projective space Pk is denoted 
simply P. 

2.1 Principalization of ideals 

Let yjo, V'l, • • • , V'p be non-zero analytic functions on a smooth complex or 
real analytic variety X. The indeterminacy points of the rational map 

ip:X—^ W 

can be eliminated as follows [I3l [2] 

Theorem 3 (Hironaka desingularization) There exists a smooth ana- 
lytic variety X and a proper analytic map tt : X — > X such that the induced 
map (f = if o 71 is analytic. 



Let Ox be the sheaf of analytic functions on X and consider the ideal sheaf 
/ C Ox generated by ipo, ipi, . . . , ipp. The inverse image ideal sheaf of / under 
the map n : X ^ X will be denoted n*I. This is the ideal sheaf generated 
by the pull-backs of local sections of I. We note that vr*/ may differ from 
the usual sheaf-theoretic pull-back, also commonly denoted by vr*/. A simple 
consequence of Theorem [3] is the following 

Corollary 1 The inverse image ideal sheaf vr*/ is principal. 

This is called the principalization of I. Indeed, as the induced map (p is 
analytic, then for every A G X there exists j, such that the functions <^j/<^j, 
i = 1,2, are analytic in a neighborhood of A. Therefore there is a 

neighborhood U of X such that (pj\Q divides in the ring of sections Ofy 
of the sheaf Ox, that is to say Ijy is generated by (pjlo. 



X 
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2.2 The Weierstrass preparation Theorem 

Definition 1 Let [a, 6] C M. A Weierstrass polynomial in a neighborhood of 
[a, b] X {0} C K X K" is an analytic function of the form 

P{u, X)=u'^ + ai{X)u'^~^a2{X) + ■■■ + ad(A) 

such that P{u,0) has exactly d zeros in [a,b] (counted with multiplicity). 

In the case a = 6 = 0, K = C, the above definition coincides with the usual 
definition of a Weierstrass polynomial in a neighborhood of the origin in C""*"^ 

m 

Theorem 4 (Weierstrass preparation theorem) Let f{u,X) be an ana- 
lytic function in a neighborhood of [a, b] x {0} C K x K" such that f{u, 0) is 
not identically zero. Then f has an unique representation f = P.h where P = 
P{u, A) is a Weierstrass polynomial in a neighborhood of [a, 6] x {0} C K x K*^, 
and h = h{u, X) is an analytic function, such that h{u, 0) 7^ 0, Vm G [a, b]. 

The proof of the above theorem is the same as in the usual case a = 6, K = C 

m 

We are also interested in the behavior of the zeros u = u{X) of the Weier- 
strass polynomial. For this reason we consider the discriminant A (A) of 
P{u, A) with respect to u. It is an analytic function in a neighborhood of the 
origin G K" which might be also identically zero. This may happen for in- 
stance if for every fixed A, such that ||A|| is sufficiently small, the polynomial 
P{u,X) has a double zero u{X), which is then analytic in A. The analytic 
function {u — m(A))^ then divides the Weierstrass polynomial P{u,X), the 
result being also a Weierstrass polynomial. These considerations generalize 
to the following 

Corollary 2 The Weierstrass polynomial P from Theorem\^has a represen- 
tation P = Pl^ .P^ . . . Pi'' , where for each i, Pi is a Weierstrass polynomial 
in a neighborhood of [a, b] x {0} C C x C", with non-vanishing identically 
discriminant Aj(A) ^ 0. 

It follows that the bifurcation locus of the zeros {u : P{u, A) = 0} of a 
Weierstrass polynomial P(., A) (and hence of an analytic function /(., A)) in 
a neighborhood of a compact interval [a, b] is a germ of an analytic set in a 
neighborhood of the origin in C". Namely, in the complex case K = C, this 
bifurcation locus B is the union of discriminant loci of Pi(., A) and resultant 
loci of pairs Pi{., X), Pj{., X). 
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In the real case K = C the real bifurcation locus B is contained in the real 
part ^{B'^) of the complex one and is also a real analytic set. The complement 
to the bifurcation locus B'^ of real zeros of a real analytic function /(., A) is 
therefore a semi- analytic set. The more general case, when f{u, 0) = on 
[a, b] will be considered in section [3l It will follow from the proof of Theorem 
[T]that the complement to the bifurcation locus B^ is then a sub- analytic set, 
see also Caubergh [3]. 

2.3 The curve selection Lemma 

Lemma 1 Let U be an open neighborhood of the origin in M" and let 

fli ■ ■ ■ 1 fki Qli ■ ■ ■ 1 Qs 

be real analytic functions on U such that the origin is in the closure of the 
semi- analytic set: 

Z := {x eU : fi{x) = ■ ■ ■ = fk{x) and gi{x) > 0,i = 1, . . . , s} 

Then there exists a real analytic curve 7 : [0, 5) ^ t/ with 7(0) = and 
7(t) EZ,yte (0,5). 

We refer to Milnor's book [18] for a proof. He does it in the algebraic category, 
but his proof works in general with minor (obvious) modifications. 

Let / = f{u, A) be a function, real analytic in a real neighborhood of 
[a, 6] X {0} C M X M" (we do not suppose that /(., 0) ^ 0). 

Definition 2 The cyclicity Cyd{{[a,b], f{.,0)), f{., X)) o//(.,0) on the in- 
terval [a, 6] C M «s the smallest integer N having the property: there exists 
£0 > and a real neighborhood V of [a, b], such that for every A G M", such 
that II A II < sq, the function f{u, A) has no more than N distinct zeros in V , 
counted without multiplicity. 

When there is no danger of confusion we shall write Cycl{[a,b], f{., X)) in- 
stead of C?/c/(([a, 6], /(., 0)), /(., A)). The number of the zeros of /(.,0) on 
the interval [a, b] C M, counted with multiplicity, is an upper bound for 
Cyd{[a,b], f{., X)), but not necessarily an exact bound. The Weierstrass 
preparation theorem. Corollary [2] and the Curve selection Lemma imply 

Theorem 5 Let f{u,X) be a real analytic function in a neighborhood of 
[a,b] X {0} C M X M" which is non-identically zero on [a, 6] x {0}. There 
exists an analytic curve [0; 5) ^ M" : £ 1-^ '^(^); -^(0) = 0, such that 

Cycl{\a, &],/(., A)) = Cycl{[a, b], /(., X{e))). 
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2.4 Cyclicity of period annuli and the Bautin ideal 

In this section the base field is K = R. 

Definition 3 Let Xx be a family of analytic real plane vector fields de- 
pending analytically on a parameter X G (M",0), and let K (Z B? he a 
compact invariant set of Xxq. We say that the pair {K,Xxg) has cyclicity 
N = Cycl{{K, Xxo), Xx) with respect to the deformation Xx, provided that 
N is the smallest integer having the property: there exists Eq > and a 
neighborhood Vk of K, such that for every X, such that \\X — Ao|| < ^o, the 
vector field Xx has no more than N limit cycles contained in Vk- If K is 
an invariant set of Xx^ (possibly non-compact), then the cyclicity of the pair 
{K,Xxo) with respect to the deformation Xx is 

Cyd{{K,Xxo),Xx) = sup{Cycl{{K, Xxg), Xx) : K (Z K,K is a compact }. 

The cyclicity Cycl{{K , Xxq), -^x) therefore the maximal number of limit 
cycles which tend to as A tends to 0. To simplify the notation, and if 
there is no danger of confusion, we shall write Cyd{K, Xx) on the place 
of Cycl{{K, Xxq), Xx). Let 11 C be an open period annulus of a plane 
analytic vector field Xq. There is a bi-analytic map identifying 11 to S*^ x A 
where A is a connected open interval. Therefore Xq has an analytic first 
integral u induced by the canonical projection 

X A ^ A : {(p,u) ^ u 

which parameterizes a cross-section of the period annulus 11. Let u i— > P{u, X) 
be the first return map and 6{u, X) = P{u, X) — u the displacement function 
of Xx- For every closed interval [a, b] C A there exists > such that the 
displacement function 6{u, X) is well defined and analytic in {{u, X) : a — eq < 
u < h + Eo, \\X\\ < Eq}. For every fixed A there is a one-to-one correspondance 
between zeros of A) and limit cycles of the vector field Xx- This allows to 
define the cyclicity Cyclijl, Xx) in terms of the cyclicity of the displacement 
function 5{u,X) on the cross section A (Definition [2]): 

Cycl{K,Xx) = Cycl{[a,hlKA)) (6) 
where K = x [a, h] (we identified 11 and 5*^ x A) and 

Cycl{n,Xx)= sup Cyc/([a,6],5(.,A)). (7) 

[a,fe]cA 

Let Mo G A and let us expand 

oo 

(5(m, A) = ^ai{X){u - uq)\ 

i=0 
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Definition 4 (Bautin ideal |21| . [20| ) We define the Bautin ideal X of 
Xx to be the ideal generated by the germs ai of ai in the local ring Oo{W^) of 
analytic germs of functions G M". 

This ideal is Noetherian and let (^i, (^2, ■ ■ ■ , <^p be a minimal system of gener- 
ators, where p = dim^I/Ail, and is the maximal ideal of the local ring 
OoiW"). Let ipi, ip2, . . . ,(pphe analytic functions representing the generators 
of the Bautin ideal in a neighborhood of the origin in M". 

Proposition 1 (Roussarie, [20]) The Bautin ideal does not depend on the 
point uq G a. For every [a,b] C A there is an open neighborhood U of 
[a, b] X {0} m M X and analytic functions hi{u, A) in U , such that 

p 

6iu,X) = J2v^WH^,^)■ (8) 

The real vector space generated by the functions hi{u,0),u G [a,b] is of di- 
mension p. 

3 Proof 

To prove Theorem [T] we may apply Theorem [5] to the displacement function 
S{u, X) ... provided that 6{u,0) is not identically zero. This is certainly not 
the case. To overcome this difficulty we principalize the Bautin ideal and 
divide the displacement map 6{u, A) by a suitable analytic function (which 
does not affect its cyclicity). 

Suppose that the cyclicity Cycl{Il, Xx) is finite. There exists an invariant 
compact subset K of Xq such that Cycl(Il, Xx) = Cycl{K, Xx). If we identify 
n to 5*^ X A, then K is identified to x [a, b] where [a, b] C A. It follows 
from ([6]), ([7]) and the definition of cyclicity, that there is an open interval a, 
[a,b] C a C A, and a convergent sequence (A'^)fc in R" which tends to the 
origin in W\ and such that for every k, 6{., A*^) has exactly Cycl(Il, Xx) zeros 
in the interval a. 

If e > is sufficiently small, then the generators ipi of the Bautin ideal, 
defined by (JHl), are analytic in the set X = {A G C" : ||A|| < £:}. According 
to Theorem [3] the rational map 

can be resolved. The projection tt : X ^ X is a proper map which implies 
that there is a convergent sequence (A'^')j in X, such that 

7r(A'=') = A'^S lim A'^' = A° G 7r-^(0) 
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and hence 

Cyd{U,Xx) = Cycl{S'x[a,b],Xx) 

= Ci/c/((5^x[a,6],X^(^o)),X^(^)) 

= Cyc/(([a,6],5(.,7r(A°)),5(.,7r(A))). 

In other words, at A° G X, the cychcity of the open period annulus is max- 
imal. Of course X^^^o) = Xq, (5(.,7r(A°)) = S{.,0) = 0. Let A be a local 

variable on X in a neighborhood of Aq, 7r(A) = A. By Corollary [T] the inverse 
image of the Bautin ideal sheaf is principal. Let o tt be a generator of the 
ideal of sections in a neighborhood of A°. By Proposition [H in a suitable 
neighborhood U^^ of [a, b] x {Aq} in M x X holds 

6{u,7TCX)) = M^)h{u,~X) (9) 

where 

ifQ = o Tc, h{u, A) = h{u, 7r(A)) 

and 

/i(M,7r(A°) = h{u, 0) ^ 0. 

We conclude that 

Cydiiia, b], 7r(A°)), 7r(A))) = Cyc/(([a, 6], 7r(A°)), 7r(A))) 

and by Theorem [5] there exists an analytic curve e t-^ A(£:), A(0) = A°, such 
that 

Cydiiia, b],hi., 7r(A°)), M-, 7r(A))) = Ci/c/(([a, 6], 7r(A°)), 7r(A(£)))). 
The curve e i-^ A(£) = 7r(A(e)), A(0) = 0, is analytic which shows finally that 

CydiU,Xx) = CydiU,Xxie)). 

Theorem [T] is proved. □ 

The proof of Theorem [2] is similar: we resolve the map (p and principalize 
the Bautin ideal. Of course we can not use the Weierstrass preparation 
theorem, neither the curve selection lemma. Let A° G 7r~^(0) be a point in a 
neighborhood of which the closed period annulus has a maximal cyclicity 

CydiU^X,) = Cyc/((n,X^(3^o)),X^(3^)). 

The displacement map takes the form Q. Suppose that A = [0, 1] param- 
eterizes a cross section to 11. Then hiu, Aq) = hiu, 0) ^ is analytic on 
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[0, 1). First of all we have to prove that the Dulac expansion of h{u, 0) is not 
identically zero. Indeed, it follows from [2ll Proposition A 2.1, p. Ill] that 
the vector field Xx possesses an analytic first integral in a neighborhood of 
n. Therefore Xx is a Hamiltonian vector field with respect to a suitable area 
form. Let 

e:£^A(£),A(0) = A° 

be any analytic curve, not contained in the zero locus of ^o{X). The displace- 
ment map of Xx restricted to the curve ^ = vr o ^ is 

S{u, nC\{e))) = c.e^{h{u, 0) + eR{u, e)),c^Q (10) 

which shows that h{u,0) = M^{u) is the Poincare-Pontryagin function as- 
sociated to the curve ^. It follows from Roussarie's theorem [21], Theorem 
C] that the Dulac expansion of h{u, 0) ^ can not be zero provided that 
0) ^ 0, and that the cyclicity of the loop 7 of Xq with respect to the 
deformation X^^^^-^-^ is bounded by the generalized multiplicity of the zero 
M = of 0). Strictly speaking, the Roussarie's Theorem is proved for the 
case k = 1 m (fTOl) . but the proof in the case k > 1 \s exactly the same. The 
fact that h{u, 0) has a Dulac expansion ^ follows also from [TOl, [9] where its 
monodromy is computed. 

We conclude that the cyclicity of the closed period annulus is bounded 
by the total number of zeros Z{M^), counted with generalized multiplicity. 
This completes the proof of Theorem On 

4 Concluding Remarks 

Theorems [T] and [2] seem to belong to the mathematical folklore. The authors 
of [5] for instance used them to compute the cyclicity of open or closed 
period annuli of particular quadratic systems with a center (but provided 
wrong references to |l5j and to the author's paper [H], see [El Remark 2.1 
and Lemma 2.1]). Particular cases of Theorems [I] and [2] were previously used 
in [H p.223-224] and [IH p.490-491]). 

As in the Introduction, let Xx and 11 be an analytic family of analytic 
vector fields and an open period annulus of the field Xxq respectively. Let 
Aix = -Mxxi^Xo) be the set of all Poincare-Pontryagin functions ^ 
associated to germs of analytic curves ^ : e ^{^), -^(0) = Aq. The set Aix 
is not always a vector space, but spans a real vector space of finite dimension 
(bounded by the number of generators of the Bautin ideal). Recall that all 
functions are defined on a suitable open interval A. We denote by Z{M^) 
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the number of the zeros of on A (counted with multiphcity) . It follows 
from the proof of Theorem [T] that 

sup Z{M^) < oo (11) 

M^eMx 

if and only if for all & M.x holds 

Z(M^) < oo. (12) 

As explained in the Introduction, the Roussarie's conjecture [201 p.23] (if it 
were true) would imply that (fTTl) . ffT2l) hold true. 

Suppose now that Xx, A G A„, is the family of polynomial vector fields 
of degree at most n and denote 

Z{n,Xx,)= sup Z{M^). 

According to Theorem [T] the number Z(n, X^q) is just the cyclicity of the 
open period annulus 11 

Z{n,Xx,) = CydiU,Xx). 

If Xxf) is a generic Hamiltonian vector field with a center, the set A4x is a 
vector space of Abelian integrals. The weakened 16th Hilbert problem, as 
stated by Arnold [T], asks to compute explicitly the number Z{n,Xxo). In 
this case the inequality ([T2l l follows from the Varchenko-Khovanskii theorem 
[23l[l7]. Suppose now that Xxq is a given plane vector field of degree n with 
a center (not necessarily Hamiltonian). The natural generalization of the 
Arnold's question is then 

Find the numbers Z{n,Xxo), Aq € A„. 

To prove the finiteness of the number 

sup Z{n,Xxo) 

then would be a generalization of the Varchenko-Khovanskii theorem. 

To answer the above questions, it is necessary to compute first the space 
of all Poincare-Pontryagin functions 

Mx = 7Wx,(Xao),A e A„. 

In the case n = 2 this is a result of Iliev [ISl Theorem 2, Theorem 3]. It is 
known that if Xq is a quadratic Hamiltonian field, then Z(2,Xo) = 2, except 
in the Hamiltonian triangle case, in which Z(2,Xo) = 3, see [TTl [5]. 
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